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ABSTRACT 



Aims. The predictions of the ellipticity of the dark matter halos from models of structure formation are notoriously 
difficult to test with observations. A direct measurement would give important constraints on the formation of galaxies, 
and its effect on the dark matter distribution in their halos. Here we show that galaxy-galaxy flexion provides a direct 
and potentially powerful method for determining the ellipticity of (an ensemble of) elliptical lenses. 
Methods. We decompose the spin-1 flexion into a radial and a tangential component. Using the ratio of tangential-to- 
radial flexion, which is independent of the radial mass profile, the mass ellipticity can be estimated. 
Results. An estimator for the ellipticity of the mass distribution is derived and tested with simulations. We show that 
the estimator is slightly biased. We quantify this bias, and provide a method to reduce it. Furthermore, a parametric 
fitting of the flexion ratio and orientation provides another estimate for the dark halo ellipticity, which is more accurate 
for individual lenses Overall, galaxy-galaxy flexion appears as a powerful tool for constraining the ellipticity of mass 
distributions. 
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1. Introduction 

The structure of cluster and galaxy halos predicted by N- 
body simulations of Cold Dark Matter show several im- 
portant features, e.g. the universal radial density profile 
(Navarro et al. 1996, 1997), and the highly non-spherical 
structure fitted well by a triaxial density profile (Jing & 
Suto 2002; Law et al. 2009). These features are related to 
the nature of dark matter as well as to the formation of 
galaxies and clusters (Kuhlen et al. 2007; Debattista et al. 
2008; Bett et al. 2010), which suggest that accurate esti- 
mates of the halo properties can provide several constraints 
on cosmology. In particular, numerical simulations with 
different assumptions predict halos with different shapes, 
e.g., simulations with non-interacting cold dark matter pre- 
dict that halos are triaxial prolate ellipsoids (Allgood et al. 
2006). Moreover, if there exists a significant difference be- 
tween the shape of a galaxy and the shape of its total mass 
distribution, this would provide additional strong evidence 
for the existence of dark matter (Suyu & Halkola 2010). In 
other words, it allows us to test how reliable the light is as 
a tracer of the matter distribution. 

Gravitational lensing provides a powerful tool for study- 
ing the mass distribution of clusters of galaxies as well as 
galaxy halos (see Bartelmann & Schneider 2001; Refregier 
2003; Schneider 2006; Munshi et al. 2008, for reviews on 
weak lensing). This is because gravitational lensing probes 
the matter distribution regardless of whether it is luminous 
or dark. The weak lensing technique has been used for clus- 
ter mass reconstructions (e.g. Clowe et al. 2006; Bradac 
et al. 2008), and for the ellipticity of the dark matter dis- 
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tribution (Cor less et al. 2009; Deb et al. 2009; Howell & 
Brainerd 2010). A method to determine galaxy halo ellip- 
ticity by stacking galaxies was proposed in Natarajan & 
Refregier (2000) and has been used to determine the ellip- 
ticities of both cluster- or galaxy-size halos (Evans & Bridle 
2009; Mandelbaum et al. 2006). A mean ellipticity of 0.46 
is found from a sample of 25 clusters using Subaru data 
(Oguri et al. 2010). 

Flexion has been recently studied as the derivative of the 
shear, and responds to small-scale variations in the gravi- 
tational potential (Goldberg & Natarajan 2002; Goldberg 
& Bacon 2005; Bacon et al. 2006). Different techniques 
were developed to measure flexion (Irwin & Shmakova 2006; 
Okura et al. 2007; Massey et al. 2007; Schneider & Er 2008). 
It was noted that flexion can contribute to studies of the 
dark matter halos in galaxies and clusters, especially for 
the detection of mass substructure (Bacon et al. 2009; Er 
et al. 2010). Flexion has been implemented on small sets of 
observational data. Okura et al. (2008) performed flexion 
measurement on data from the Subaru telescope to detect 
substructure in Abell 1689. Leonard et al. (2007) analyzed 
images taken by the HST Advanced Camera for Surveys, 
from which a preliminary galaxy-galaxy flexion signal has 
been detected. Hawken & Bridle (2009) studied galaxy halo 
ellipticity using flexion. It is found that the constrains from 
flexion are comparable to, or even tighter than those from 
shear. Moreover, since flexion drops off faster than the shear 
as one goes away from the center of the lens, multiple de- 
flections from the three-dimensional mass distribution be- 
tween us and distant source galaxies (Howell & Brainerd 
2010) may not be significant for flexion. 
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In this paper, we present a new approach to estimate 
dark matter halo ellipticity using galaxy-galaxy lensing flex- 
ion. The spin-1 flexion, a vector, is decomposed into its 
two components, the radial and tangential flexion. The 
tangential-to-radial flexion ratio yields an estimate of dark 
halo ellipticity. We need to assume that the center of the 
galaxy or the cluster be known. Whereas for galaxies this 
may be a lesser problem, the mass centroid of clusters is 
sometimes difficult to determine. We study the effect of 
noise coming from an intrinsic flexion. In Sect. 3, we per- 
form a series of numerical tests of our estimate and study 
how our results are biased by the distribution of background 
sources, the centroid offsets and intrinsic flexion. We discuss 
our results in Sect. 4. 



2. Basic formalism 

The full formalism described here can be found in Bacon 
et al. (2006); Schneider & Er (2008). Weak lensing shear 
and flexion are conveniently described using a complex no- 
tation. We adopt the thin lens approximation, assuming 
that the lensing mass distribution is projected onto the 
lens plane. The dimensionless projected mass density can 
be written as n(9) = E(0)/£ cr , where 8 is the angular po- 
sition, £(<?) is the projected mass density, and S cr is the 
critical surface mass density. 

The first-order image distortion by gravitational lensing 
is described by the shear 7, which transforms a round source 
into an elliptical image. The higher-order effect, which is 
called flexion, is described by two parameters. The spin-1 
flexion is the complex derivative of k 

^ „ 9k , dn .„ . 

and the spin-3 flexion is the complex derivative of 7, 

Q = V c7 • (2) 



2.1. Radial and tangential flexion 

The spin-1 flexion is a vector-like quantity, the gradient 
of the surface mass density. Therefore, the spin-1 flexion is 
directed towards the center of the lens in the case of an axi- 
symmetric mass distribution. However, if the lens deviates 
from axial symmetry, e.g., due to elliptical halos or mass 
substructure (Bacon et al. 2009; Hawken & Bridle 2009), 
the flexion vector will have a different direction. 

In galaxy-galaxy lensing, shear can be decomposed into 
tangential and cross components. Analogously, we decom- 
pose the spin-1 flexion into a radial and tangential flexion 
component. They are defined as 



-T-r 
-T ■ 4> 



(3) 
(4) 



where the spin-1 flexion vector is given by T = {T\,T<2) = 
Vk in Cartesian coordinates, f and are the unit direction 
vectors in polar coordinates. For the spin-3 flexion, there 
is no such a clear intuitive picture of its two components; 
thus, we only consider the spin-1 flexion in this paper. 

As mentioned before, the spin-1 flexion vector is not 
directed towards the center of an elliptical lens (Fig. 1). 
Therefore, the tangential flexion no longer vanishes and can 
be used to estimate the ellipticity of the mass distribution. 



2.2. Elliptical mass distributions 

We now assume that the isodensity contours of the mass 
distribution are ellipses with ellipticity e, or equivalently, 
axis ratio (1 — e)/(l + e), and orientation O of the major 
axis. In this case, the surface mass density can be written 
as 



k{0) =K 6U 



0o) , sin 2 (0-0 o ) 



' cos^ 



(1 + c) 2 



[l-ef 



(5) 



where the function K{p) describes the radial density profile. 
The derivatives of k with respect to radial and azimuthal 
coordinates are given as 




cos 2 (0 - O ) sin 2 (0 - O ) 



(1 + e) 2 (1-e) 2 
ldn _ K' 2e(l - e 2 )- 2 sin2(0- O ) _ 



cos 2 (0— 0o ) . sin 2 (0— 4> ) 
2 1 



= ~Tr , (6) 
T T ■ (7) 



Both of these derivates have the same radial profile, given 
by K' . We now define the flexion ratio as 



•Ft 



Tr 



(8) 



which is the tangent of the angle between the direction to 
the mass center and the direction of the flexion vector. For 
an elliptical mass distribution, this becomes 



r(0) 



|2esin2(0-0 o )| 
1 - 2e cos 2(0 - 0o ) + e 2 



(9) 



Thus, the flexion ratio r(0) depends only on the elliptic- 
ity e and the orientation 0o of the mass distribution, not 
on its radial profile. In particular, it is independent of the 
lens strength, and thus of the source and lens redshift. In 
principle, from the measurement of the flexion ratio at two 
polar positions 0, one can determine both e and 0o, though 
due to noise, this determination will have large uncertainty. 
Alternatively, one can consider the flexion ratio averaged 
over all polar angles, 



(r) 



1 

2^ 



d0 







2e sin 2(0 — 0o) 



1 - 2e cos 2(0- 0o ) + e 2 



2, He 

= -In- 

tt 1 — e 



(10) 



which depends solely on the ellipticity e. Due to its simplic- 
ity, this mean flexion ratio can be easily measured from a 
flexion field in a given aperture. 

However, we have to assume to know the center of 
the lens for calculating the two flexion components. For 
a galaxy, the center is assumed to be the bright center of 
galaxy. For clusters, the location of the BCG not necessarily 
coincides with the mass center. Moreover, 'intrinsic flexion' 
also introduces extra noise (see below). 

For a given value of e, the flexion ratio is bounded. One 
can see from (9) that 



r{4>) < r max = 



2e 



1-e 2 ' 



(11) 



Arbitrary large ellipticities, which would allow r to become 
very large, are implausible and most likely do not exist. We 
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Fig. 1. Spin-1 flexion vector field for an elliptical isothermal 
density distribution with e = 0.3. T only points towards the 
center when the background galaxy is located on the major or 
minor axis of the halo. Note that the length of the vectors are 
logarithmically scaled, for better visibility 

can take limits on e from other observations, e.g., Parker 
et al. (2007) compared galaxy-galaxy lensing along the ma- 
jor and minor axes, and found that the axis ratio b/a of 
galaxy halos lies between 0.5 and 0.8 (e G [0.11,0.33]). 
Moreover, Oguri et al. (2010) analyzed 25 clusters and 
found (1 — b/a) — 0.46, corresponding to e = 0.32. Here 
we use the conservative assumption that e < 0.8, putting 
an upper bound to the flexion ratio of r < 4.5. We will 
employ this limit to remove excessively large flexion ratios 
which are due to noise. 



3. Numerical test with NIE toy model 

In this section we describe some simulations which we have 
performed in order to test the behavior of the estimator 
given in the previous section. For this, we model the surface 
mass density profile by a non-singular isothermal elliptical 
profile (NIE), described by 



K(p) 



(12) 



where 6>e is the Einstein radius, describing the strength of 
the lens, and C is the core radius; for C = 0, this special- 
izes to the singular isothermal elliptical profile (SIE). In 
our simulations, we take 9e = 6" and 8 C = 2". Galaxies are 
randomly distributed within an 1' x 1' 'source plane' behind 
the lens. Resulting images are discarded if they are located 
closer to the halo center than 6" (the strong lensing regime) 
or at distances \0\ > 30", where the flexion signal will be 
very small. Moreover, very large flexions cannot be mea- 
sured (Schneider & Er 2008), thus images with \T\ > 0.5 
are discarded as well. 

The flexion ratio is calculated according to Eq.(8) for 
each point. A mean flexion ratio is obtained by 



1 N 



(13) 



?=i 



where N is the number of images for each realization. Then 
the mass ellipticity is estimated, according to Eq.(10), as 

e= L —^r. (14) 
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Fig. 2. Comparison of the ellipticity estimator (14) with the 
input ellipticity (solid line). The estimated ellipticity for 20 re- 
alizations each for 20 input values of e are shown by pluses, for 
NIE models with 9 C = 2" 



3.1. Noise- free case 

Due to the non-linearity, the estimator (14) is expected to 
be biased. In order to test this, we generated mock data sets 
with different e = 0.03z, i = 1,2,. ..20. For each ellipticity, 
we used 20 realizations, with a density of flexion points 
of 80arcmin -2 . The filtering described above leads to 63 
flexion data for each realization on average. 

For each realization, we calculated f from Eq.(13) and 
estimate e using Eq.(14). In Fig. 2, we show e vs. the input 
value of e. The solid line shows the identity, the plus points 
are the estimates from the 20 realizations. In fact, the esti- 
mates e closely trace the input value. The variance between 
different realizations becomes larger with increasing e. The 
mean over the 20 realizations is shown by the plus points 
in Fig. 3, showing that e very slightly underestimates the 
input value. 



3.2. Intrinsic flexion 



gTTf/2 + I ■ 



Our main source of noise is intrinsic flexion, meaning that 
the sources can have non-vanishing third-order brightness 
moments. Depending on the method of flexion measure- 
ment, the intrinsic noise may be different (Goldberg & 
Leonard 2007). Since flexion has the dimension of an inverse 
length, the intrinsic flexion is also inversely proportional to 
the image size. Therefore, the distribution of intrinsic flex- 
ion depends on the survey and is difficult to obtain from 
real measurements on current data. Effects of a point spread 
function (PSF) need to be considered, which is far more dif- 
ficult for flexion than for shear measurements. In particular, 
an anisotropic PSF may affect the direction of the flexion 
vector. 

We use a simple model to generate intrinsic noise for our 
simulated data, by setting F ohs = T\ + n/i + i(J r 2 + n/2). 
The components of the intrinsic flexion n/i,n/2 are drawn 
from a Gaussian distribution, with each component being 
characterized by ajr\ — ajr 2 — 0.03 arcsec -1 . 

In Fig. 3 we compare the estimates e with the input val- 
ues (solid line). The plus points show the mean over 20 
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realizations without noise (Fig. 2), and the dashed line dis- 
plays estimates with intrinsic noise included. Not surpris- 
ingly, our estimates from noisy simulations are larger than 
the input values (intrinsic noise will cause the estimate to 
deviate from zero even for a perfectly symmetric mass dis- 
tribution). This is particularly significant for small e. As an 
additional test, we reject flexion data which have r > 4.5, 
and perform our estimate again. The result is shown by 
crosses in Fig. 3. One can see that the bias is slightly re- 
duced in this case, in particular for larger e, but still e over- 
estimates the true ellipticity. The amplitude of this bias 
depends of course on the intrinsic flexion distribution. 

Another way to estimate the ellipticity of the mass dis- 
tribution is to fit the flexion ratio with the model (9), with 
the two free parameters e and </>o. The result for two real- 
izations are presented in Fig. 4. The points are calculated 
from simulated data of an NIE with e = 0.3 (left) and 
e = 0.6 (right). The dotted lines are the fitting result with 
Eq.(9) whereas the solid curves show Eq.(9) with the input 
ellipticity. The fits almost perfectly agree with the input 
model. For the left panel, the fitting yields e = 0.31 and 
</>o = 0.1, whereas e = 0.36. In the right panel, e = 0.62 and 
4>o = 0.006 from fitting, whereas e = 0.67. In both panels 
of Fig. 4, there are several points with significant deviations 
from the fitting curves which is the reason for e to be larger 
than the input value, whereas these points do not strongly 
affect the fitting result. 

Obviously, the fitting method yields a more accurate 
estimate of the mass ellipticity than the estimator e, and at 
the same time also estimates the orientation. It is therefore 
the preferred method for individual lenses. In contrast, the 
estimator e can be applied also to an ensemble of lenses, by 
superposing their respective flexion values. In this way, e 
estimates the weighted mean ellipticity of the ensemble of 
lenses. Note that in contrast to galaxy-galaxy lensing with 
the shear method, the individual lenses do not have to be 
aligned before the averaging; hence, no assumption about 
the relative orientation of mass and light needs to be made. 



3.3. Centroid offset test 

Since the ellipticity represents the deviation from symmetry 
of a mass distribution, a centroid offset will affect the deter- 
mination of the ellipticity, leading to an additional bias. We 
simulated this effect, by calculating the radial and tangen- 
tial flexion components with respect to 6q + 56, where 0q is 
the true center of the lens, and 56 is the offset. For simplic- 
ity, we choose two sets of 59, one along the major axis, the 
other along the minor axis of the mass distribution. To iso- 
late this effect, we generate 400 data sets without intrinsic 
noise, for two NIE with e = 0.3, 0.6. The results are shown 
in Fig. 5. A centroid offset indeed biases the estimate of 
the ellipticity, but the magnitude of this effect is relatively 
small, as long as the offset is considerably smaller than the 
Einstein radius of the lens. For the lens with higher ellip- 
ticity (bottom panel), an offset along the minor axis has 
a larger impact on the resulting ellipticity than one along 
the major axis. We conclude that the effect of misidentify- 
ing the lens centroid is of little concern for galaxy lenses, 
but may be more important for galaxy clusters where the 
centroid is often ill-defined from observations. 
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Fig. 5. Bias due to the centroid offset. The solid line is our 
estimate (400 sets of data) without centroid offset. The plus 
(cross) points are estimate for the centroid offset alone the major 
(minor) axis of the elliptical halo. The top (bottom) panel is for 
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4. Conclusions and Outlook 

In this paper, we have studied galaxy-galaxy flexion with 
elliptical mass distributions. We derived the ratio r of radial 
and tangential flexion and found that it is independent of 
the radial mass profile of the lens, and of the source and lens 
rcdshift. The flexion ratio depends solely on the ellipticity 
and orientation of the mass distribution. We defined an es- 
timator e for the ellipticity in terms of the flexion ratio, 
and tested its performance with simple simulations. This 
estimator does not rely on knowing the orientation of the 
distribution, i.e., assuming that it follows the orientation 
of the light; hence, it can be used to statistically superpose 
several lenses and thereby obtain their average ellipticity. 
The independence of this estimator from an assumed ori- 
entation is quite different from methods to constrain lens 
cllipticities from shear (Braincrd & Wright 2000; Hoekstra 
et al. 2004), for which the orientation needs to be known. 
In the appendix we show that an analogous method to that 
presented here does not work for shear, but is unique to the 
spin-1 nature of the first flexion. 

For an individual lens with several flexion measure- 
ments, the ellipticity and orientation can be derived by 
fitting the measured flexion ratios with these two parame- 
ters. The accuracy of the resulting estimate is substantially 
better than that of e. 

Our simulations showed that e is biased. A first bias is 
due to the non-linearity of the relation between estimated 
ellipticity and mean flexion ratio. The size of this bias de- 
pends on the number density of sources for which flexion 
can be measured, and disappears in the limit of very high 
source density (or large number of lenses for which the mean 
ellipticity is explored). Intrinsic flexion causes a more sig- 
nificant bias, in particular for small values of e. A misiden- 
tification of the mass centroid is another source of bias; 
however, if the centroid offset is much smaller than the 
Einstein radius of the lens, this bias is rather small. 

There will be several complications in a real analysis. 
First, the flexion T can not be directly measured, but only 
reduced flexion (Schneider & Er 2008). This only requires 
a small modification, though: if k is constant on ellipses, so 



will be ln(l — n). Since 



V c ln(l - k) = 



Gi - ggj 
1 - 99* 



(15) 



where G\ is the spin-1 reduced flexion, G\ = V*<? = 
(J 7 + gT*)/{l — k), and g = 7/(1 — n) the reduced shear, 
we can use the phase of (15) instead of the flexion ratio. 
Furthermore, the impact of intrinsic flexion is uncertain, 
reflecting our lack of knowledge about its magnitude which 
is ill-determined from current observations. A further care- 
ful study of intrinsic flexion will be of interest before the 
flexion method is applied to real data. 

In addition, as we have shown, a low number den- 
sity of background sources also introduce bias. Thus a 
deep survey with more background sources is clearly 
favoured. Current surveys such as the Hubble Space 
Telescope (HST), the Canada- France-Hawaii Telescope 
Legacy Survey (CFHTLS) or Subaru telescope might be 
able to constrain halo ellipticity using flexion, if the flexion 
can be measured with sufficient accuracy. The James Webb 
Space Telescope will almost certainly allow very accurate 
measurements of the ellipticity of galaxy- and group-sized 
lenses. 
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Appendix A: The shear profile in SIE halo 

In this appendix, we present the tangential and cross shear 
profile for an SIE halo model and show that one cannot 
apply the ratio of cross to tangential shear to estimate the 
ellipticity of the lens. 

To calculate the shear, we calculate the deflection po- 
tential using the Poisson equation 



ld_ 

6d0 



+ 



1 d 2 ip 



2k . 



(A.l) 
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The convergence can be written in polar coordinates 
0e(1 - c) 1 



K{6,<t>) 



29 



vT 2 



cos + sin 



(A.2) 



where / is the axs ratio. The solution of the Poisson equa- 
tion can be found in Kormann et al. (1994), but here we 
write it in our notation as 



0E(l-e)fl 



(A.3) 



arcsin (/ cos 0) cos + arcsh ( — sin J sin 



where /~ = v 1 — f 2 - Then we rewrite the potential in 
Cartesian coordinates and calculate the shear by 71 = 
(-0,11 — V',22)/2 and 71 = -0,12- We thus obtain the two 
shear component for the SIE halo, 



71 = 

72 = 



E (l-6)(0?-flj) 

20 2 d 
E (1 - e)M2 



6 2 d 



(A.4) 
(A.5) 



9 2 . According to the definition 7 t 



where d = \J f 2 6\ 
—71 cos 20 — 72 sin 20, 7 X = 71 sin 20 — 72 cos 20, the 
tangential-to-cross shear for SIE halo reads 



7t = 



c(l-e) . 
2d ' 



0. 



(A.6) 
(A.7) 
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